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ABSTRACT 



Experimental data for the charging of an air receiver is 
presented and interpreted in detail. The data indicates a 
substantial departure from the adiabatic behavior. The ex- 
perimental results are used to evaluate existing closed form 
expressions for the thermodynamic state of a gas in a re- 
ceiver. A method for experimentally determining the convec- 
tive heat transfer coefficient is developed, evaluated and 
used in conjunction with these expressions. 

The experimental work was performed from March 1969 
through May 1969 at the Naval Postgraduate School, Monterey, 



California. 
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I . 



INTRODUCTION 



Systems using the charging of a receiver are often em- 
ployed today. A common practice when analyzing a receiver 
being charged is to assume that the process is adiabatic. 

This simplifies the governing equations for the system and 
leads to a fairly straightforward expression for the state of 
the gas in the vessel. During the course of the charging 
process however, a substantial temperature difference between 
the walls of the receiver and the gas in the receiver may 
develop. Often the tank's thermal capacitance is of suffi- 
cient magnitude to permit the extraction of large amounts of 
energy from the gas with only a small change in the tempera- 
ture at the outside surface of the tank walls. This obscures 
the apparent effect of heat transfer, for the absence of a 
noticeable temperature change at the outside of the receiver 
walls may lend credence to the adiabatic assumption whereas 
actually large amounts of heat are being transferred from the 
gas to the receiver walls. Errors introduced by such an 
assumption may lead to serious problems. For example, if 
heat transfer is neglected, calculations for determining the 
amount of gas needed to charge a vessel av lead to a pre- 
dicted value lower than that actually required. 

In order to aid in the analysis of this problem Reynolds 
[Ref. 1] developed a theory which includes the effects of 
heat transfer for determining the thermodv '^mic state of gas 
in a receiver during the charging process. He developed four 
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closed form solutions for the various magnitudes of the system 
parameters. Using his criteria to determine the proper closed 
form solution, many charging processes may be analyzed with- 
out excessive difficulty. The accuracy of these solutions de- 
pends to a great extent on the accuracy with which the average 
convective heat transfer coefficients can be estimated for the 
charging process. 

In his experimental investigation of the blowdown process 
[Ref. 2] Reynolds tested his theory using an h (average con- 
vective heat transfer coefficient for the entire blowdown 
process) determined by taking the value of his heat transfer 
parameter NTU (number of thermal units, see Appendix A) which 
when used in his theoretical equation of state gave the best 
fit to the experimental data. These h values were on the 
order of those predicted on the basis of an assumption of 
steady state turbulent free convection inside the receiver, 
i . e . , 

^ = . 13 (GRPR) 1//3 

Therefore Reynolds suggested the use of an h based on this 
assumption in his solution to the charging process. 

The problem still remains, however, that unless h is de- 
termined by an independent means it is difficult to come to 
any quantitative conclusions as to the accuracy of Reynolds' 
theory . 

The purpose of this study has been to derive a method for 
experimentally determining the average convective heat trans- 
fer coefficient between the gas and the receiver, and to 
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apply these results to three of Reynolds' closed form solu- 
tions in order to evaluate their accuracy. The three closed 
form solutions evaluated are for conditions where the heat 
transfer from the receiver walls to the surrounding medium 
can be neglected. These three cases are the ones most com- 
monly encountered. The fourth of these closed form solutions 
deals with the case where the walls of the receiver are very 
thin and the heat transfer from the receiver walls to the am- 
bient medium must be taken into account. This solution was 
not evaluated due to the difficulties in experimentally deter- 
mining the average convective heat transfer coefficient be- 
tween the walls of the receiver and the surrounding medium. 
Also presented in this paper are representative values for 
the various heat transfer parameters determined in this 
undertaking . 

In the following sections of this thesis the theoretical 
and analytical methods and results are discussed. Following 
a description of the charging system, the theoretical model 
due to Reynolds is reviewed and the methods for the experi- 
mental evaluation of the heat transfer coefficient are devel- 
oped. The experimental results are then discussed in the 
light of their relevance to Reynolds' model and their useful- 
ness in engineering problems. A final section is included in 
order to summarize this work and indicate some areas for 
future investigation. 
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II. OBJECTIVES 



Experiments on the charging of a gas receiver were con- 
ducted at varying flow rates and heat transfer environments, 
in order to accomplish the following objectives: 

(1) Evaluate the closed form solutions developed by 
Reynolds to approximate the thermodynamic state 
of a gas during charging. 

(2) Develop, use, and evaluate an expression for ex- 
perimentally determining the average instantaneous 
heat transfer coefficient in a gas receiver being 
charged. 

Ill . DESCRIPTION OF EXPERIMENTAL APPARATUS 
A. GENERAL DESCRIPTION 

The test apparatus (Figures 1 through 4) consisted of an 
aluminum cylindrical gas receiver with a volume of 1.03 cubic 
feet. The complete physical dimensions of the test apparatus 
are listed in Table 1. The vessel was designed according to 
the ASME Boiler and Pressure Vessel Code for a working pres- 
sure of 200 psig. Its top was removable to allow access to 
the interior of the tank, using an O-Ring flange arrangement 
to insure airtightness when the top was in place. The re- 
ceiver was fitted with a flow metering device mounted in the 
center of the tank top. The other fittings consisted of a 
mounting for a thermocouple probe, a pressure transducer 
mounting, and a bleed valve arrangement. 
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This test receiver was mounted so as to be immersed in a 



larger tank which was either open to the atmosphere or filled 
with ice water depending on the particular run. Air for 
charging the test tank was supplied from two 117 cubic foot 
air vessels charged by an air-cooled compressor. The supply 
pressure to the flow metering device was varied from 180 psig 
to 10 psig. 

B . FLOW METERING 

The mass flow rate for each run was determined by maintain- 
ing a critical flow through one of two nozzles or through an 
orifice. These small diameter (D) devices were machined from 
stainless steel in accordance with a paper by Grace and Lapple 
[Ref. 3] and the discharge coefficients (C D ) given in their 
paper were used in these calculations. The diameter and dis- 
charge coefficients are listed in Table 1. The orifice or 
nozzle was mounted on the tank top in a flange type arrange- 
ment that was sealed by 0-Rings and connected to a 1 inch line 
that lead to the supply tanks. The stagnation pressure (Pq) 
was read by a local pressure tap and gage arrangement mounted 
just upstream of the flow metering device. The stagnation 
temperature (Tq) was likewise found by using a thermocouple 
mounted in the flow just upstream of the metering device. It 
should be noted here that for velocities of the magnitude oc- 
curring in these tests, it was safe to assume that the differ- 
ence between stagnation and local conditions upstream of the 
orifice was negligible. Thus, assuming adiabatic flow through 
the orifice or nozzle, the temperature and pressure values of 
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P Q and T q indicated by the apparatus just described can be 
used in the critical flow equation. 



w = .532 C 



D 






The air to the flow metering device was controlled by a cruick 
acting Jamesbury ball valve. 



C. TEMPERATURE MEASUREMENT 

The test vessel was fitted with a thermocouple probe con- 
sisting of four 40-gage copper constantan thermocouples. The 
thermocouple housings were designed so as to measure as 
closely as possible the local temperature at four equally 
spaced heights in the receiver (Fig. 1) . The thermocouples 
were arranged so that the four could be read in series, and, 
at the same time, any one of the four could be read independ- 
ently. The series signal was recorded on channel one of a 
two channel Moseley strip chart recorder, while the signal of 
the thermocouple being read singly was recorded on a contin- 
uous Brown recorder. According to the manufacturer's speci- 
fications thermocouples of this type are accurate to within 
±1.5°R for the temperature ranges found in these tests. 

Since the thermocouple probe readings were only used qualita- 
tively, no further calibration was made. As was mentioned 
above, the stagnation temperature of the inlet air was read 
from a single thermocouple. The signal of this 40-gage copper 
constantan thermocouple was monitored on a Leeds and Northrup 
recorder. 
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D. PRESSURE MEASUREMENTS 



The pressure in the test receiver was obtained by use of a 
Daystrom Wiancko pressure transducer attached to the receiver 
whose signal was recorded on channel two of the Moseley strip 
chart recorder. The input voltage (21 volts) to the trans- 
ducer was supplied by a Philco power supply and the transducer 
output signal was adjusted so as to read 10 osi/in on the 
Moseley recorder. The pressure transducer was calibrated 
originally using a test gage and later using a dead-weight 
tester. When connected to the recorder the pressure trans- 
ducer registered pressures within ± 1 psi of the actual pres- 
sure throughout the range of these tests. 

IV. EXPERIMENTAL PROCEDURE 

The experimental program was divided into the following 
three cases : 

1. In the first group of runs the receiver was charged 
12 times under conditions corresponding to the ad- 
iabatic case in Reynolds' theory. The runs con- 
sisted of three charges at each of four possible 
supply tank pressures; 190, 180, 170 and 160 osig. 

The 3/8 inch knife-edge orifice wa mounted for 
these tests and the total elapsed time for a run 
was approximately one second. The large tank en- 
circling the test receiver was open to the atmos- 
phere. Thus the rapid chargings gave NTU values 
which dictated the use of Reynolds ' adiabatic 
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model for theoretically determining the thermodvnamic 
state of the gas in the receiver throughout the run 
(See Section V-A-3 ) . 

2. The second group of runs corresponded to the isother- 
mal case described in Reynolds' theory. The 36 runs 
consisted of 3 runs at each of 12 possible supply 
tank pressures ranging from 180 psig to 10 psig, in 
various intervals. The 1/32 inch nozzle was mounted 
in the flow metering apparatus and the run time var- 
ied from 150 seconds to 75 seconds. Once again the 
large tank surrounding the test tank was open to the 
atmosphere. The resulting very small mass flow rates 
were introduced to correspond to NTU values in the 
region of 7, a lower limit for Reynolds' criteria. 

3. The last set of runs was made with the test receiver 
completely immersed in ice water contained by the 
large tank. The first 9 runs were made with the 
1/32 inch nozzle in place, three each at supply pres- 
sures of 180, 160 and 140 psig. The last 18 runs of 
this case were made with the 1/8 inch nozzle mounted 
in the tank top, three each at supply pressures of 
180, 150, 120, 90, 60 and 30 psig. The ice bath was 
to maintain the receiver walls at a constant tempera- 
ture and thus present the gas with an isothermal sink 
to correspond to the one described in the analytical 
model . 
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The values of the supply pressure listed above represent 
the nominal charging pressure. The actual staanation pressure 
just upstream of the flow metering device was recorded for 
each run. 

In making a series of runs, the supply tanks were first 
charged to approximately 190 psig. The system was then iso- 
lated from the compressor and the test receiver was allowed 
to reach equilibrium with its surrounding medium. The quick 
acting valve between the supply tanks and the flow metering 
device was then opened. It was kept open until the pressure 
in the test tank reached a value insufficient to maintain a 
critical pressure ratio across the flow meter. The valve was 
then closed and a run terminated. For the succeeding run the 
supply tank pressure was adjusted to the desired value and 
the test receiver was vented to the atmosphere. The next run 
did not commence until the temperature of the gas in the re- 
ceiver (as indicated by the thermocouple probe) reached equi- 
librium with that of the tank walls and thus the surrounding 
medium. Having made one pass through the supply pressures 
indicated for a particular case, the supply tanks were re- 
charged and the procedure repeated for a second and third 
time . 

Throughout the charging process the pressure and tempera- 
ture history of the gas in the tank was recorded on a two 
channel strip chart recorder. 
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V. THEORETICAL CONSIDERATIONS 



A. THEORY DEVELOPED BY REYNOLDS 1 
1 . Introduction 

In an analysis of the charging process one object is 
to obtain the time dependent thermodynamic state of the aas 
in the receiver. A common practice is to express the tempera- 
ture of the gas (the dependent variable) as a function of the 
mass (the independent variable) . For a constant volume re- 
ceiver the perfect gas law may then be used to express the 
pressure as a function of the mass. Knowing the mass flow 
rate the thermodynamic state of the gas as a function of time 
can be described. 

With this in mind Reynolds developed the following 
model for a receiver being charged: 




The thermal capacitance of the receiver walls and any other 
internal structure is lumped into a single capacitance 



x The nomenclature used in this section is taken from Ref- 
erence 1, see Appendix A for listing. 
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represented by C c and equal to the product of the mass of the 
capacitance and its specific heat (C c = M c c ) . The heat trans- 
fer resistance between the receiver walls and the surrounding 
environment is represented by and the heat transfer resist- 
ance between the capacitant material and the gas in the re- 
ceiver is represented by R^. 

2 . General Differential Equation for Charging 

Reynolds assumes that the thermal resistance as well 
as the thermal capacitance are invariant with time. He also 
assumes perfect mixing of the injected gas and the gas in the 
receiver, and that the walls of the receiver are at a uniform 
temperature throughout. Using these assumptions he extends 
his model to include all the significant enerqv terms as 
follows : 




Combining an energy balance for the above model and the heat 

transfer rate equations, 

T -T T-T 

c «> c 

q «> " R q i " R. 

oo 1 

the following non-dimensional general differential equation 
for the charging is obtained (see Appendix B) : 



2i 



NTU+NTU 



U * M * 35^- + w * t2 “* + NTU + M * §£ + - c;*"“ “‘’Sf 



. * , NTU NTU 

+ [W * dM^ + c-T(NTU + NT Ua5 )w* + - - g -g — °° ] T* 

0 0 

, * NTU+NTU NTU NTU 

- kT i w * as* ~ ' c 0 * ~ ** kT ! - c 0 T ~ = 0 • 

For the case of constant mass flow rate (w* = 1) this equation 
reduces to: 



d 2 T* rlT* 

M* — -f- + [C, +C 0 M* ] 2i_ + C^T* + C. = 0 
dM* 2 1 2 dM* 3 4 



( 1 ) 



where , C 2 , C^, and are constants comprised of the sys- 
tem parameters NTU, NTU and C^*; defined as, 



(hA) 



NTU = 



NTU = 



R . c w rt w c ' 

1 v 0 o v 



(hA) 



00 R -o c v W 0 *0^ 



and 

Me C 

C * = C --P = — 

0 M.c M.c 
0 v 0 v 

The parameter NTU is used by Reynolds to represent the conduct- 
ance (reciprocal of resistance) of the receiver as well as to 
give a measure of the rate of the process. The Cq* parameter 
is the ratio of the capacitance of the receiver to the initial 
capacitance of the gas. 

Reynolds solves equation (1) analytically, however 
the resulting series solution is seen to be impractical for 
many engineering applications. Fortunately in many cases the 
values of the system parameters are such that a simpler closed 
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form solution may be obtained. These solutions can be obtain- 
ed by simplification of the series solution or by returning 
to the general differential equation for charging. In his re- 
port Reynolds nonetheless chooses to develop each closed form 
solution from a simplified form of his original model. A 
brief description of these and their solutions are presented 
in the following four sections. A complete derivation of the 
general differential equation for charging as well as deriva- 
tions of the closed form solutions are presented in Appendix 
B. 

3 . Adiabatic Charging 

As has been mentioned, a common method of analysis of 
a system being charged is to assume the process is adiabatic. 
This assumption of no heat transfer to or from the gas may be 
found to represent the case of a high mass flow rate charging 
(Wq large) process fairly accurately. This may be thought of 
as due to the fact that the gas has not had time to transfer 
a substantial amount of heat to the walls before the charging 
process is terminated. Even if the high mass flow rate proc- 
ess is fairly lengthy, the effect of the heat transfer is 
small compared to the effect of introducing large quantities 
of mass during the charge. Such an assumption may also be 
justified for the case where the gas has been insulated from 
the thermal capacitance of the receiver walls (h - 0) . This 
case corresponds to the situation where NTU is verv small, 
and the following model applies : 



23 



dU = d (Mu) 



/////////A 

dM 

/7777777V/ 

/ M,T 

/ 77777777777777777 



TV 

'/ 
'/ 



The solutions for the dimensionless temperature and 
are , 



T* 



kT* 



kT*-l 

"M* 



pressure 

( 2 ) 



and 



P* = kT*(M*-l) + 1 



( 3 ) 



4 . Isothermal Charging 

A receiver with a large thermal capacitance (C 
large) and either charged very slowly (Wq small) or having a 
high hA value, will exhibit an isothermal behaviour. The 
normal temperature increase in the gas due to compression is 
not observed because the energy is extracted from the gas bv 
the thermal capacitance of the receiver before it can become 
significant. As described above, a system of this type will 
have large values for the system parameters NTU and Cq*, and 
may be represented by the following model : 
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The temperature ratio for charging is by definition, 

T* = 1 



and therefore, 



P* = M* 



5 . Charging at Constant Mass Flow with Heat Transfer 
to an Isothermal Sink 

In some systems the thermal capacitance of the wall 
is much larger than the thermal capacitance of the aas (Cq* 
very large) . For this case it can reasonably be assumed that 
the energy extracted from the gas is not great enough to sig- 
nificantly affect the temperature of the receiver walls (T c 
is a constant) . Thus the analytical model is modified to 
describe the heat transfer between the gas a'nd an isothermal 
sink as follows: 




The solutions for this case are, 



T* = 



kT*+NTUT* - (kT?-l-NTU+NTUT* ) M* (1+NTU ) * 

1 . c 1 c 



1+NTU 



(4) 



and 



P* = 



(kT*+NTUT* ) M* - (kT*-l-NTU+NTUT*)fl* 
1 c 1 C__ _ 

1+NTU 



-NTU- 



• (5) 
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6 . Charging at Constant Mass Flow with Inside Resistance 
Negligible 

A system having a high conductance can be assumed to 
have negligible inside heat transfer resistance (NTU high) . 
This implies that the temperature of the capacitance is the 
same as that of the gas and the heat transfer is between the 
capacitance and the surrounding environment. That is distin- 
guished from the isothermal case in that here the thermal 
capacitance is finite but not large (very thin walled cylinder 
for example) . Since systems of this nature are not too com- 
mon, Reynolds feels that the greatest value of this solution 
is that it supplies information for determining the effect of 
capacitance on the heat transfer in a charging process. The 
model for this system is: 




The solutions for this case are 




1+NTU 



00 



+ kT* + NTU T* 



00 00 



( 6 ) 



1+NTU 



00 
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and 



1 




1+NTU 




00 



P* = T* (C*+l) 



(1+NTU ) T*-kT*-NTU T* 
' 00 2 _ 00 00 



(7) 



B. CRITERIA FOR APPLICATION OF CLOSED FORM SOLUTIONS 

In the previous four sections, closed form solutions were 
developed for the state of the gas in a receiver being charged. 
These solutions were developed by simolifving the original 
model with assumptions as to the magnitude and imoortance of 
various terms in the general solution. For example, in the 
adiabatic charging case the parameter NTU is assumed to be 
very small, thus the mass flow rate is large compared to the 
heat transfer coefficient and a solution based on no heat 
transfer is formulated. Similar types of assumptions are 
made for the other three cases. Reynolds then produces quan- 
titative criteria for the use of such assumptions. These 
criteria are based on a maximum deviation of 5% in T* at the 
value of M* = 9. They are determined by comparing the case 
in question with the case which would give the maximum depar- 
ture from this behavior. Thus the adiabatic case, (no heat 
transfer to the walls) is compared to the case where the ca- 
pacitance of the receiver walls is infinite, the isothermal 
sink case (large amounts of heat transfer to the walls). 

When the percentage T* difference between the two solutions 
is plotted against NTU, the value of the two solutions agree 
within 5% for values of NTU less than .25. In a similar 



fashion criteria are developed for all the closed form solu- 
tions. Some solutions have more than one possible criteria 
for their use, however, only the criteria for the three cases 
tested here which are met by our exoerimental system are 
listed below. As mentioned before, the three closed form 
solutions evaluated here are the most useful and the criteria 
for their application are those most commonly found. 

Adiabatic: 0 < NTU < 0.25, all 

Isothermal Sink: 0.25 < NTU <7 > 40 

Isothermal: NTU >7 C* > 40 

C. DETERMINATION OF CONVECTIVE HEAT TRANSFER COEFFICIENT 

In order to test Reynolds' theory on the state of a gas 
in a receiver being charged, a technique was devised to obtain 
a value for the convective heat transfer coefficient in the 
receiver in which the tests were run. To this end the follow- 
ing model of the experimental system was constructed: 



H.m 



l 
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where , 



m = the mass flow rate into the receiver 
q = the heat transfer to the receiver walls 
d(Mu)/dt = the rate of change of the internal energy 
of the gas in the receiver with resoect 
to time 

= the stagnation enthalpy of the entering gas 
T = the temperature of the gas in the receiver 

P = the pressure of the gas in the receiver 

M = the mass of the gas in the receiver 

Several points should be noted here. First it was assumed 

that the temperature throughout the receiver could be repre- 
sented by a single value T (See discussion of imperfect mixing 
section V-E) . The same was true for the pressure term P. It 
should also be noted that this development does not include 
the case where there is heat transfer between the surrounding 
medium and the walls of the receiver. 

An energy balance on the receiver yields, 

d (Mu) 



mH i - q = 



dt 



The kinetic energy of the gas in the receiver was neglected 
for it can be shown that for the flow rates in this exoeri- 
ment, the velocity of the gas in the receiver was verv small 
when compared to the internal energy. Using the notation f 
to represent the derivative of a function f with respect to 
time, the above equation can be written as. 



- q = 



mu + Mu 
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Assuming the gas in the receiver is thermallv and calorif- 
ically perfect, 

q = -m(c T-c T.) - Me T (8) 

M v p 1 v 

where is the inlet stagnation temperature of the gas. By 

definition of the convective heat transfer coefficient, 

hA(T-T ) = q - -mc T - m(c T-c T.) 
w M V V d l 

In the same manner as T represents the average gas temperature, 
T w is used to represent the average temperature of the re- 
ceiver walls. Thus the h defined here represents an average 
instantaneous heat transfer coefficient for all points on the 
inside of the receiver walls. The symbol A represents the 
total inside area of the receiver walls. By use of the per- 
fect gas law it can be seen that, 

P _ m T 

P M T 



for a constant volume system. Therefore the above equation 
can be written as , 

-Mc T(2-)+mc T. 
v, _ V P pi 

h A(T-T ) 

w 

or, using the perfect gas law again. 



h = 



-c P (^) +riic T . 
v R pi 

[ (£Y) _t ] 

Al W V 



( 9 ) 



The next step was to determine an expression for T w if it 
was not maintained a constant as in the isothermal sink case. 
In doing so the difficult task of determining experimentally 
an accurate average wall temperature was avoided. Since, for 
simplicity, no heat transfer from the outside of the receiver 
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walls to the surrounding medium was assumed, an energy balance 
on the walls yields , 

q = isr(M t u ) = M CT 
at w w w w 

where C is the specific heat of the walls and is assumed con- 
stant. Combining this with equation (8), 

q = -m(c T-c T.) - Me T = M CT 
M v p 1 v w w 

is obtained. Integrating this equation over time from t = 0 

to t = t gives , 

M CT dt = j!f _ )r[-m(c T-c T ) -Me T]dt 
J t=0ww J t= 0 V PI V 



or 



M C(T -T ) = ft n -4r-(Mc T) dt + „ me T.dt 

w w w Q J t=0 dt v J t=0 p l 



-t=t 



In these tests T was a constant, therefore, 

M C(T -T ) = -Me T + M n c T„ + (M-Mjc T . 
w w Wq v 0 v 0 0 p l 

where subscript 0 denotes the condition at time t = 0. 
Solving for gives 



T = 
w 



-Me T+M n c T_+Mc T . — M n C T. 
V OvO pi Opi 

M C 
w 



+ T 



w. 



Again using the perfect gas law to eliminate the term T, the 



final form of the expression for T is obtained, 

w 



T = T + 

W Wq 



M n c \ .. 
! 0 o) 



M C 
\ w 



[ ir T i- TJ ' 

Hq ± 



c v V 
Irm c 

\ w 



[p-p 0 i 
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Equations (9) and (10) furnish a means for estimating the 
film coefficient h in terms of the tank pressure and gas mass 
and their rates. The massive experimental simplicity thus 
introduced is the major justification for the acceptance of 
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inaccuracies stemming from the definition of the heat transfer 
in terms of spacial averages of gas and wail properties. 

D. DISCUSSION OF ASSUMPTION OF NO EXTERNAL HEAT TRANSFER 

In the derivation of the expression for h it was assumed 
that there was no heat transfer from the outside of the re- 
ceiver walls to the surrounding medium. This simplifies the 
energy balance, for if there is no heat transfer from the walls 
to the outside medium, all the energy leaving the gas must be 
stored in the receiver walls. This assumption is not onlv 
valid for this experimental system, but for many pressure ves- 
sels used in charging processes. A vessel capable of with- 
standing a substantial pressure is likely to be constructed 
with materials and dimensions that give it a fairly large 
thermal capacitance when compared to that of the entering 
gas. In the test receiver, for example, the thermal capac- 
itance of the receiver was approximately 8.7 Btu/°R while 
even at its highest mass the gas capacitance was only of the 
order of .08 Btu/°R. Thus the gas temperature drop due to 
heat transfer to the receiver walls would have to be extremely 
large before any significant change in the wall temperature 
was observed. 

This assumption of no external heat transfer may also be 
justified by considering the length of time it would take the 
temperature increase due to heat transfer at the interior sur- 
face to reach the outside walls. One method of analyzing this 
transient conduction heat transfer problem would be to consider 
the receiver walls a serai-infinite solid with heat transfer 
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due to convection at the surface. Thus by determining the 
time (t) required for an observable temperature change to 
occur in the semi-infinite solid at a distance (x) from the 
surface equal to the receiver wall thickness of 3/8 inch, some 
idea of the temperature change at the outside surface of the 
test tank will be obtained. According to Carslaw and Jaeger 
[Ref. 4], the solution for this type of oroblem is given by 
the equation, 



v 

V 



erfc 



( x \ hx+h 2 at 
2/at ) 



erfc (-^ — 
\2 /at 




where 



a = the thermal diffusivity of the material 
h = h/k = the convective heat transfer coefficient/ 
the thermal conductivity 

v = the difference between the wall temperature at 
time t and a distance x from the surface and the 
initial wall temperature 

V = the difference between the gas temperature and the 
initial wall temperature. 

Using the highest value of V recorded throughout our runs 
(100°R) and a corresponding value for h we find that v is 
still only approximately 1.5°R even after 2 minutes have 
elapsed. Thus for our runs the assumption of no external heat 
transfer during a charge is justified. This will be the case 
for all receivers of thermal capacitance sufficient to store 
the thermal energy, transferred from the gas, with negligible 
increase in temperature. 
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E. DISCUSSION OF THE EFFECTS OF IMPERFECT MIXING 

In the derivation for determining h, and in the closed 
form solutions derived by Reynolds, it was assumed that the 
effects of imperfect mixing could be neglected. In other 
words, it was assumed that for analytical purposes the gas 
can be considered perfectly mixed so that no temperature or 
pressure gradients exist. As was seen in an unpublished ex- 
periment on the heat transfer in a closed container (with 
similar dimensions to the test receiver) after gas injection 
[Ref. 5], this is not actually the case. Temperature gra- 
dients do exist and this leads to uncertainty in evaluating 
the heat transfer data correctly. In order to determine the 
effective temperature potential for heat transfer, use must 
be made of some sort of average gas temperature so as to keep 
the complexity of the experimental equipment and the data re- 
duction process reasonable. If this average gas temperature 
is dependent upon the degree of mixing, then the assumption 
that the gas is perfectly mixed may lead to erroneous 
conclusions . 

In his paper [Ref. 1], Reynolds showed through the use of 
a simple comparison between mixed and ummixed gas systems that 
the pressure is independent of the degree of mixing in an ad- 
iabatic receiver regardless of its magnitude. The magnitude 
of heat transfer in these tests was never very large and, 
since the average gas temperature used in the derivation for 
determining h is based directly on the perfect gas law, the 
temperature thus calculated is fairly independent of the de- 
gree of mixing. 
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Reynolds also showed that the temoerature averaged with 
respect to mass in such an adiabatic receiver was independent 
of the degree of mixing. Realizing that this average temoer- 
ature would be very difficult to find experimentally, he went 
on to compare the temperature averaged with respect to mass 
and that averaged with respect to volume and found that the 
discrepancy was so small "that use of a volume average tem- 
perature in experimental investigation is entirely satisfac- 
2 

tory." Therefore, the average temperature measured by the 
four thermocouples connected in series at equal volumes 
vertically in the tank could be used as the temperature refer- 
red to as the gas temperature in the closed form solutions. 

VI . SUMMARY OF EXPERIMENTAL RESULTS 
A. FORM OF RESULTS 

The data from a particular run consisted of a pressure and 
temperature history for the test receiver during charging, and 
values for the inlet stagnation pressure and temperature. 

From this data, using the critical flow equation, the mass 
flow rate for each run was calculated and thus the mass of the 
gas could be determined for any instant of time. Using this 
mass and the pressure history of the run, the instantaneous 
h values were determined by equations (9) and (10). The val- 
ues of NTU were then calculated and used to determine which 
class of closed form solution the criteria would designate 



'Reference 1, p. 



85. 
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for this run. The values of P* = P/Pq as determined bv the 
appropriate closed form solution were then compared with the 
experimental values. The average percentage deviation be- 
tween these two values and the average value for NTU through- 
out a run were then tabulated. Plots were made comparing the 
P* values as predicted by Reynolds and those observed experi- 
mentally for representative runs. These results are in the 
form of P* versus M* = M/Mq because the pressure readings 
were slightly more accurate due to the uncertainties involved 
in determining the temperature of a gas in motion. 

The values of h used in determining NTU have been repre- 
sented by plots of the Nusselt number (NU = hL/k^) versus 
the product of the Grashof and Prandtl (GRPR = I^p^gB^AT/ 
y | (CpVi/k) f) numbers for various runs. The subscript f here 
indicates that the term was evaluated at the film temperature, 
which for this case was the mean temperature between that of 
the gas and the wall. The characteristic length L of the 
system was taken to be the receiver height for these 
calculations . 



B. HIGH MASS FLOW RATE RUNS 

A series of runs was conducted with the test receiver ex- 
posed to the ambient atmosphere at high mass flow rates so as 
to correspond to the adiabatic closed form solution developed 
by Reynolds. The dimensionless capacitance term 
C* = C c /MqC v , used in the analytical criteria was approx- 
imately 657 for all runs. Based on the criteria NTU < .25 
for any value of C* , all of these runs were found to correspond 
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to the so called "adiabatic" case. Values of 1.4 for the 
ratio of specific heats of the gas (air in this case) and 1.00 
for the non-dimensional inlet stagnation term T* , (the inlet 
stagnation temperature of the gas, / the original tempera- 
ture of the gas, Tq) were used in calculating the theoretical 
P* values from equation (3). Figure 5 is a graphical compar- 
ison of the actual and theoretical P* values for representa- 
tive runs of this type. NU versus GRPR plots for a few 
typical runs of this nature are seen in figures 11, 12, 13 
and 14 , the points having been plotted at 1/10 second inter- 
vals throughout the runs. Table 2 summarizes the results of 
this series of tests. 

C. LOW MASS FLOW RATE RUNS 

A second group of runs was conducted at low mass flow 
rates so as to correspond to the isothermal closed form solu- 
tion proposed by Reynolds. Once again, the test receiver was 
exposed to the surrounding atmospheric conditions. Reynolds' 
criteria for an isothermal charge solution to hold is NTU > 7 
and C* > 40. The value of C* for this group of runs was 
again approximately 657 and the NTU values calculated varied 
from 1.8 to 14.6. Thus these runs straddle the cutoff value 
of NTU = 7 and can be used to determine the validity of this 
figure. The solution for an isothermal charge is simply 
P* = m* , for T is a constant and therefore T* = 1. The val- 
ues of P* predicted by Reynolds and those observed experi- 
mentally are compared graphically in figure 6 for represent- 
ative runs of this type. The h values found in a few 
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typical runs are presented by the NU versus GRPR graphs of 
figure 15, 16, 1/ and 18. The points in these figures repre- 
sent 15 second intervals throughout the runs. Table 3 sum- 
marizes the results of this series of tests. 

D. ISOTHERMAL SINK RUNS 

The final group of runs was made, with the test receiver 
immersed in an ice bath, at intermediate mass flow rates so 
as to correspond to Reynolds' isothermal sink model. The 
criteria for this case is simply C* > 40. However, in order 
to avoid crossing into the adiabatic or isothermal solution 
regions, the value of NTU must be greater than .25 and less 
than 7 respectively. It should be understood that the ice 
bath was not necessary in order to meet Reynolds' criteria: 
with or without it the value of C* (approximately 629 for this 
case) was large enough to dictate the isothermal sink solution. 
The ice bath was employed to ensure that the temperature of 
the receiver walls was a constant and that a significant tem- 
perature potential existed between the receiver walls and the 
gas. Values of 1.4 for the ratio of the specific heats of the 
gas and 1.00 for the dimensionless inlet temperature, Tj* , as 
well as for the dimensionless wall temperature, T*, (the re- 
ceiver wall temperature T / the original temperature of the 
gas Tq) were used in calculating the theoretical P* values 
from equation (5) . Figure 7 is a graphical comparison of the 
actual and theoretical P* values for representative runs of 
this type. NU versus GRPR plots for a few typical runs of 
this nature are seen in figure 19, 20 and 21. The points 
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were determined at 15 second intervals for figure 19, and 1 
second intervals for the other two figures. Table 4 sum- 
marizes the results of this series of tests. 

VII. DISCUSSION OF RESULTS 

In order to systematically analyze the test results, the 
following section is divided into three separate sections. 

The first section is devoted to an analysis of the equation 
derived for the convective heat transfer coefficient. The 
results obtained when comparing the experimentally observed 
values of P* to those predicted by Reynolds are then discus- 
sed. The last section deals with an analysis of the convec- 
tive heat transfer coefficients obtained with emphasis on the 
identification of general trends that might be investigated 
in a more comprehensive study on this topic alone. 
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A. ANALYSIS OF THE EXPRESSION FOR DETERMINING THE CONVECTIVE 
HEAT TRANSFER COEFFICIENT EXPERIMENTALLY 
In an earlier section the exoression 



h = 



* V* • 

-c P (=•) +mc T . 
v R pi 

A[(— )-T 1 

U MR ; "vv 



( 9 ) 



was derived with the value for T , when it was not considered 

w 



a constant, given by equation (10). 



T = T + 
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As was mentioned before, these expressions were puroosely 
developed to obtain a value for h from experimentally deter- 
mined data. Due to the difficulty of measuring the effective 
average temperature of the gas, this term was eliminated in 
the above expressions. The equation for h is thus a function 
of P and M, the other terms being system oarameters which were 
generally held constant and known to a high degree of accuracy 
throughout a run. 

First, considering the expression for the wall temperature, 
equation (10) , it is seen that even when the temperature of 
the receiver was not held constant the temperature of the walls 
changed very little during a run. This can be seen by noting 
the magnitude of the terms in parentheses as well as the fact 
that the pressure and mass of the gas are never more than one 
order of magnitude larger than their original conditions (Mq, 
Pq) . An investigation of the values calculated for T^ shows 
that this was indeed the case for these tests. 
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Next the expression for the heat transfer coefficient it- 
self will be considered. Since the derivation of this expres- 
sion was centered about the fact that, 

h = A(T?T ) 
w 

it can be seen that the numerator of expression (9) is an ex- 
pression for the heat transfer q and the denominator is the 
product of the inside area of the receiver walls and the tem- 
perature potential between the walls and the gas. Considering 
first the denominator, it can be seen that when the tempera- 
ture of the gas and that of the walls are the same, this term 
is zero. Of course in this situation, the heat transfer (the 
numerator of the expression) is also zero and thus as would 
be expected the convective heat transfer coefficient is un- 
defined. Now, examining the case where there is a finite but 
small difference between the wall and gas temperatures, as 
might be the case in a low mass flow rate charging process, 
it is seen that the equation for h becomes very sensitive to 
errors in the temperature of the gas or of the receiver walls. 
It has already been observed that the expression for the tem- 
perature of the walls is insensitive to small errors in the 
values of the system variables, the gas pressure and mass. 
Unfortunately this is not the case for the perfect gas rela- 
tion used to express the gas temperature. Even though an 
error in T may seem negligible when considered on a percentaae 
basis, it is not when considering small differences between T 

and T . For example, an error of 1°R in the temperature of 
w 

the gas will cut the value for the denominator of h in half 
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if the temperature difference xs only 2°R. A temperature error 
of this magnitude is definitely withxn the experimental accu- 
racy of these tests, therefore one must exercise caution when 
applying this equation to runs in which the temperature poten- 
tial between the receiver walls and the gas is small. 

The numerator of equation (9) can be analyzed in much the 
same manner. As was the case with the denominator, it can be 
shown (by examining the normal magnitude of the various terms) 
that except for the case of small temperature differences this 
expression is capable of absorbing small errors in P and M 
without producing a significant error in h. 

Combining these results, it is concluded that equation (9) 
is insensitive to small error in P and M when the temperature 
difference between the gas and the tank walls is at least of 
the order of 10°R or greater. Below this value the calculated 
values of the heat transfer coefficient may be subject to 
substantial error. 

The effect this conclusion has on the confidence in the 
results of these tests is discussed in the following sections. 
The actual value of h is used only in the eauation for the 
value of P* in the isothermal sink solution. The runs that 
correspond to this model have substantial temperature poten- 
tials throughout, therefore the values of P* can be calculated 
with confidence. In the other two cases considered, the val- 
ues of P* are not a function of the convective heat transfer 
coefficient, but rather the value of h is used in determining 
the range of NTU in order to determine the oartic alar solution 
to be employed. The values of NTU for these runs are dictated 
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to a greater extent by the mass flow rate than by the convec- 
tive heat transfer coefficient. For this reason, although the 
values of NTU may not be exact, they can certainly be used to 
evaluate the relationship between NTU and the deviation of the 
theoretical P* values from those actually observed in either 
of these cases. 

The situation is less clear cut when it comes to inter- 
preting the convective heat transfer coefficient data. The 
sensitivity of the expression for h at low values of the tem- 
perature potential places the numerical values for some runs 
in doubt. The low mass flow rate runs in particular are 
greatly affected by this uncertainty. For these runs the tem- 
perature potential was only on the order of 10°R for the ini- 
tial runs and decreased as the charging pressure was progres- 
sively lowered for each successive set of runs. This does not 
mean that this data is not meaningful. There is no reason to 
believe that the trends indicated by such data should not be 
correct; however, it must be realized that the numerical re- 
sults are subject to a large range of experimental uncertainty. 

B. COMPARISON OF REYNOLDS' MODELS WITH THE EXPERIMENTAL 

RESULTS 

1 . High Mass Flow Rate Runs 

As was mentioned earlier, all runs of this nature were 

meant to correspond to Reynolds' adiabatic solution criteria. 

The value for C* was 657 in all cases and the calculated val- 
o 

ues of NTU were well below .25 in each run. Reynolds would 
then recommend that the solution to this type of charging be 



43 



approximated by the assumption that there was negligible heat 
transfer to the receiver walls. When a comparison of the val- 
ues of the dimensionless pressure as predicted by this solution 
and those observed experimentally was made the deviation be- 
tween the two values was of the order of 7.5% of the actual 
value (see Table 2) , with the theoretical value being higher 
than the actual. 

The discrepancy between these two figures may be at least 
partially accounted for by noting that the assumption that 
there is no heat transfer is not exact. Since the effect of 
heat transfer is to lower the change in internal energy of the 
gas and thus diminish the increase in the gas temperature, a 
lower gas temperature results with a subsequent reduction in 
pressure. With this in mind one expects the theoretical P* 
value to be higher than that actually observed, as is the case 
for this data. 

The parameter NTU is a measure of the thermal conductance 
(h) relative to the rate of gas flow (Wq) . The value of NTU 
should then be related to the exactness of the assumption of 
no heat transfer. A connection between NTU and the deviation 
of the actual and theoretical P* values can be noted in a 
single run. NTU is directly proportional to h for any partic- 
ular run, for Wq is a constant. An increase in h for runs of 
this sort is accompanied by an increase in the deviation (Fig. 
5) . Although the quantitative relationship between h and Wq 
is not clear from the data for different runs, the increase in 
h and Wq is such that their ratio, and hence NTU, increases 
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with increasing flow. The expected increase in the departure 
of theory from experiment with increasing NTU is observed in 
the comparisons (Fig. 8) . 

Reynolds proposed, on theoretical grounds, that for NTU 
less than .25 the deviation between actual values of P* and 
T* and those predicted by the adiabatic solution would be less 
than 5% for values of M* as high as 9 . In the experimental 
data, this figure seemed valid at the lowest values of NTU 
(-.08) recorded and appeared slightly optimistic at higher 
values (10% deviation at NTU =.12). In any case, the accu- 
racy of these tests was insufficient to discriminate a 5% 
deviation and the criteria developed by Reynolds is suffi- 
cient for engineering calculations. 

In conclusion, it has been seen that the deviation be- 
tween the actual thermodynamic state of a gas and that pre- 
dicted by the adiabatic solution for charging is a function 
of the actual magnitude of the heat transfer in the system. 

The actual temperature and pressure were always found to be 
lower than that predicted by the adiabatic solution for the 
receiver. The parameter NTU proposed by Reynolds seems to be 
a good measure of whether or not an adiabatic solution should 

be attempted. His criteria seems satisfactory for all but 

« 

the most exacting engineering calculations as well. There is 
also evidence to the fact that the values of the convective 
heat transfer coefficient and the mass flow rate are not inde- 
pendent in charging processes of this nature. 
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2. Low Mass Flow Rate Runs 



As was discussed in an earlier section, the runs of tnis 
type were performed in order to be used in evaluating Reynolds' 
isothermal model. The value for Cq in all cases was 657 and 
the values of NTU ranged from 1.8 to 14.6. For values of NTU 
greater than 7, Reynolds recommends the use of the isothermal 
closed form solution for charging processes. When the values 
of the dimensionless pressure, P*, as predicted by the isother- 
mal solution and those observed experimentally were compared 
the deviation between the two values was somewhere between .5% 
and 4.0% of the actual value (see Table 3). In all cases this 
deviation proved to be negative; that is, the actual pressure 
and temperature were higher than those predicted in an iso- 
thermal solution to the charging process. 

The deviation in the isothermal values of P* and the 
actual values for these runs were small and can be accounted 
for by investigating the assumptions that lead to the iso- 
thermal charging model. In this model the thermal capac- 
itance of the receiver was assumed infinite while the internal 
heat transfer resistance was taken to be zero. Therefore, any 
tendency for the incoming gas to heat up was countered by an 
immediate transfer of heat to the receiver walls. This con- 
dition was approached by charging so slowly that the gas and 
the walls, initially at the same temperature, remained at that 
state. This model can not be exact since the interfaces be- 
tween the gas and the walls have a finite heat transfer resis- 
tance, thus the resulting discrepancies between the actual and 
theoretical values for P*. 
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It seems reasonable to predict in view of the preceding 
discussion that an increase in the heat transfer resistance 
(a decrease in h) would lead to an increased deviation in the 
isothermal and actual values for P*. Although the values for 
h were subject to error in these runs (see section VII-A) , a 
definite downward trend throughout a single run was noted 
(Fig. 15, 16, 17 and 18) . As predicted, this decrease in h 
was accompanied by an increase in the deviation of the P* 
values (Fig. 6) . 

An examination of the effect of NTU on the deviation of 
the isothermal and actual thermodynamic state of the gas was 
then made. It was noted that the values of h for different 
runs remained constant within experimental scatter regardless 
of the value of Wg. Thus, unlike the adiabatic case where 
the effects of changing Wg seemed to be offset by correspond- 
ing changes in h, the value of NTU for this type of charging 
was inversely proportional to the value of w Q . Once again 
NTU proved to be a good measure of the deviation, for as Wg 
decreased the percentage deviation decreased (Fig. 9) . 

The criteria proposed by Reynolds suggests that for values 
of NTU larger than 7 an isothermal solution for a charging 
process will give less than 5% deviation from the actual ther- 
modynamic state of the gas . The data from these tests indi- 
cates that for values of NTU as low as 2 the deviation is 
still less than 5%, therefore it seems that Reynolds' criteria 
is slightly conservative. Of course, it must be pointed out 
that the value of for these runs was 657 and this is well 
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above the value of C* > 40 listed along with the NTU > 7 cri- 
teria. Nonetheless, it can safely be asserted that for values 
of above 40 and values of NTU of 7 or greater, the isother- 
mal solution to a charging process will give good accuracy. 

In conclusion, it was observed that the actual temperature 
and pressure in a receiver being charged are consistently 
higher than those predicted by an isothermal solution. The 
parameter proposed by Reynolds once again appeared to be a good 
indicator as to whether or not an assumption of isothermal 
charging can be applied to a particular process. Reynolds' 
criteria also proved to be very satisfactory for engineering 
purposes. Finally, evidence was observed that the value of 
the convective heat transfer coefficient and the mass flow 
rate can be considered independent in a charging process of 
this nature. 

3 . Isothermal Sink Runs 

The last group of runs to be evaluated wore those made 
with the test receiver immersed in an ice bath. These runs 
were performed at intermediate values of NTU so as to corre- 
spond to Reynolds ' isothermal sink closed form solution. The 
only criterion specified for the use of this method is that 
the dimensionless capacitance term be greater than 40. In all 
the runs of this nature the value for was 626, therefore 
this criteria was easily satisfied. Of course, Reynolds im- 
plies that the value of NTU dictating an isothermal sink solu- 
tion (to approximate the actual conditions) be be .ween .25 and 
7, for if not the regions encountered prescribe either the 
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adiabatic or the isothermal solution. When the values of the 



dimensionless pressure as predicted by the isothermal sink 
solution and those observed experimentally were compared, the 
deviation between the two values was somewhere between 8% and 
14.5% of the experimental value. 

It should be noted here that the isothermal sink solution 
is an intermediate case between the two extremes of an adia- 
batic solution and the isothermal solution. Therefore, the 
thermodynamic state of the gas must be expressed as a func- 
tion of the heat transfer coefficient and the mass flow rate 
in the closed form solution (see equations (4) and (5)). 

As before, physical interpretations were sought for the 
deviation between the isothermal sink solution for the thermo- 
dynamic state of the gas and that actually observed. The 
assumption on which this model was built was that the temper- 
ature of the thermal capacitance remained a constant even 
after absorbing the energy released by the gas in the receiver. 
If the comparatively small amounts of heat transferred to the 
receiver walls were uniformly distributed throughout the very 
large thermal capacitance, this assumption might be quite 
accurate. However, as was seen by an investigation using the 
heat conduction equation for a semi-infinite solid (see sec- 
tion V-D) most of the internal energy (as indicated by the 
temperature distribution) never gets any farther than the 
first few tenths of an inch in a charaing orocess of this na- 
ture. For this reason the temperature at the wall surface 
does not remain constant but rather increases slightly even 
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though the effective thermal capacitance of the receiver is 
still extremely large. Bearing this in mind, the actual tem- 
perature potential between the inside surface of the receiver 
walls and the gas is overestimated by the isothermal sink 
solution. Therefore, the amount of heat transferred out of 
the gas and the accompanying temperature and pressure drops 
are slightly exaggerated in this closed form solution (Fig. 

7) . 

To begin an evaluation of the parameters that affect the 
deviation between theory and experiment, single runs, in 
which the mass flow rates were constant, were examined. For 
the low mass flow runs there seemed to be only a very gradual 
increase in the value of h as the run proceeds. At higher 
values of the flow rate the increases in h during a run were 
more substantial as were their magnitudes (Figs. 19-21). 

Using both the low mass flow runs (almost constant h) and the 
high flow rate runs (increasing h) it was possible to examine 
the effect of h on the deviation between an actual run and 
the thermodynamic state predicted for this run by the iso- 
thermal sink solution. In analyzing the data (Figs. 7 and 
10) the conclusion was reached that the deviation was appar- 
ently not closely dependent on the value of h. In both the 
low flow rate runs and the high flow rate runs the deviation 
followed the same pattern throughout a run and likewise there 
was no discernible connection between the magnitude of the 
deviation and the values for h. In any case it is probable 
that the effect, if any, of the value of h on the deviation 



10 



of experimental and isothermal sink values is masked by the 
use of NTU = hA/WQC v as a parameter in equations (4) and (5). 

For the two extreme solutions (adiabatic and isothermal) 
presented in the previous sections, their accuracy increased 
as the appropriate extreme condition was approached. This 
might lead one to believe that the intermediate ranges of NTU 
for the isothermal sink solution now being considered would 
most closely correlate to the actual thermodynamic state of 
the gas during charging. Since the data does not cover the 
entire range of NTU between .25 and 7, no conclusive statement 
can be made in this regard. Figure 10, however, supports the 
notion that as NTU approaches an intermediate value of say 3.5 
the deviation decreases. One thing can be asserted, and that 
is that, due to the connection between w Q and h in this type 
of charging, the values of NTU do not vary greatly. Thus, 
regardless of the connection between NTU and the deviation, 
it seems certain that the deviation between the isothermal 
sink solution and that actually observed is approximately 
constant for a particular system as long as the flow rates 
are of the same order of magnitude. 

In conclusion, it should be noted that the isothermal 
sink closed form solution is a function of the convective 
heat transfer coefficient and the mass flow rate, and thus 
differs from the previous two closed form solutions. The 
deviation between the thermodynamic state of the gas as pre- 
dicted by the analytical solution and that observed experi- 
mentally was between 8% and 15% of the experimental value. 

The closed form solution values were consistently lower than 
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those found m the tests. This discrepancy can be explained 
by the fact that the predicted value of the temperature of 
the thermal capacitance in Reynolds' solution is lower than 
that actually existing. The heat transfer and temperature 
drop are subsequently overestimated in this solution. For 
runs of this nature it was also observed that the values of 
h increase throughout a run and that not only the macrnitude 
of this increase but the values of h are increased when the 
mass flow rate is increased. There seemed to be little or no 
dependence of the accuracy of the isothermal sink assumption 
on the values for h when Wq is a constant. No conclusive 
statement can be made as to the connection between NTU and 
the deviation of the isothermal sink solution and the actual 
values , although there is evidence to support a prediction 
that this solution is most accurate at intermediate values of 
NTU. Due to the apparent connection between Wq and h the 
values of NTU do not change substantially for runs in this 
region nor do the percentage deviation in the values for the 
thermodynamic state of the gas. Lastly, for this particular 
experiment the discrepancy between the solution as predicted 
by an isothermal sink solution and that observed experimen- 
tally, was slightly higher than that predicted by Reynolds 
even though the value of C* was well above the criteria value 
of 40. 
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C. DISCUSSION OF THE EXPERIMENTALLY OBTAINED HEAT TRANSFER 

COEFFICIENTS 

Having empirically obtained values for the convective heat 
transfer coefficient throughout the various runs, an attempt 
was made to obtain at least a qualitative explanation of the 
behavior of these values. The charging situation was sub- 
stantially different from those usually described in convection 
theories. The pressure of the system is usually not varied 
in such theories, nor is the direction and magnitude of the 
flow velocity a function of variables other than position. 

Since any attempt to arrive at a theoretical expression for 
the convective heat transfer coefficient in a receiver being 
charged would be completely beyond the scope of this paper, 
no attempt was made. Some insight was gained, however, in 
considering the possible non-dimensional parameters that 
might affect an expression for the convective heat transfer 
in a charging situation. 

From the general equations of motion and energy, an 
average convective heat transfer coefficient evaluated over 
the total surface area of a body (this eliminates any spatial 
dependence) can be expressed as a function of the dimension- 
less parameters as listed below: 

NU = hL/k = f (PR , ER, GR, RE) 

An investigation of these parameters was then made in order 
to come to a better understanding of the causes for the behav- 
ior of h in the experimental data. First to be considered 
was the Prandtl number. Defined as 

PR = yc /k 
P 
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this dimensionless number can be considered constant and of 
the order of 1 for solutions involvinq perfect gases. Having 
concluded that the effects of the Prandtl number need not be 
considered, the Eckert number was then considered. This non- 
dimensional parameter expresses the ratio of the inertia terms 
to the enthalpy terms and is defined as 

ER = U 2 /c AT 
P 

where U is the characteristic velocity of the flow and c AT 
is the change in enthalpy across the boundary layer. In the 
charging processes carried out for this paper the average gas 
velocities were small as compared to the enthalpy, therefore 
the effect of the Eckert number is considered negligible. 

As has been pointed out, a common practice is to take the 
Grashof number as the governing parameter for the heat trans- 
fer processes of the kind under consideration here. Defined 
as , 

GR = gL 3 BATp 2 /y 2 

this term is the ratio of the buoyancy to the viscous effects. 
The buoyancy of the fluid aids the heat transfer by supporting 
the motion of the fluid at the surface of the receiver walls. 
Noting this, the Grashof number must be considered significant 
in determining the convective heat transfer coefficient in a 
charging process of the sort presented here. Since consider- 
able motion is also present in the tank due to the stirring 
action of the injection, the Reynolds number, defined as 

RE = Up£/y 
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must also be considered. In the above expression, U is again 
a characteristic velocity of the flow and & a characteristic 
length for the flow. This dimensionless number can be seen 
to be a measure of the relative importance of the inertia and 
viscous effects and often occurs in the analysis of problems 
in forced convection. 

Investigating the variables comprising these non-dimen- 
sional parameters leads to some understanding as to their 
magnitude and behavior for charging processes. For small 
temperature changes the density can be considered proportional 
to the pressure. If the temperature changes are significant 
the density is both proportional to the pressure and inversely 
proportional to the temperature. The viscosity is almost 
independent of pressure and increases slightly with tempera- 
ture. For situations of an isothermal nature, it is then ex- 
pected that the variables in the Reynolds number are the 
velocity, the density (which is proportional to the pressure) 
and perhaps, the characteristic length of the fluid flow. If 
the temperature changes throughout a process are large, once 
again the Reynolds number is proportional to the velocity, the 
density (now proportional to the pressure and inversely pro- 
portional to the temperature), the characteristic length, and 
also inversely proportional to the viscosity (which is propor- 
tional to the temperature) . Now applying the same sort of 
reasoning to the Grashof number, the effect of y as a variable 
in a moderate temperature change system can again be eliminated 
along with the gravitational constant, g. The density behaves 
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in the manner described in the Reynolds number discussion. 
The coefficient of volumetric expansion, $, is defined as 
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-I(iP) 

p ' 



P 



and varies as the inverse of the temperature for a perfect 
gas. The term AT in the Grashof number is the temperature 
change across the boundary layer and L is the characteristic 
length of the system, a constant. 

In preparation for applying what was learned about the 
factors affecting the heat transfer coefficient to the behav- 
ior observed in the experimental runs, the velocity variations 
in the test receiver were examined. When the inlet valve for 
the receiver was first opened a pressure ratio (the receiver 
back pressure/ the stagnation pressure upstream of the flow 
meter) considerably smaller than the critical pressure ratio 
existed across the flow metering device. Sonic flow was 
established at the throat of the inlet to the tank. Therefore, 
fluid entering at a pressure corresponding to the sonic condi- 
tions was forced to adjust to the lower back pressure in the 
tank. Regardless of the form of this adjustment it must have 
been coupled to a substantial increase in the local velocity. 

As the back pressure in the receiver increased, the pressure 
adjustment became less severe and thus the local velocity de- 
creased throughout a charging process. Here the term "local" 
applies to the conditions in the tank near the point of 
injection . 

For the low mass flow rate runs the temperatu. e potential 

between the gas and tne walls, AT, was very email. For this 
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reason it was felt that the Grashof number was not as signif- 
icant as the Reynolds number throughout these runs. Having 
seen that during a charge the velocity near the injector de- 
creased, it was thought that perhaps the characteristic length 
over which the incoming jet acts was also decreased. There- 
fore, although the density increased due to the increase in 
pressure, it seemed reasonable to assume that the Reynolds 
number decreased slightly throughout a run. A droo in the 
Nusselt number, and thus the convective heat transfer coeffi- 
cient, was predicted. In other words, forced convection gives 
way to free convection as the run slowly proceeds. Figures 
15, 16, 17 and 18 seem to support a conclusion of this sort. 
This notion also is supported in a paper by Ulrich et. al. 

[Ref. 6] in which the mode of heat transfer following gas 
injection was observed to proceed from forced to free convec- 
tion. For the high mass flow rate and isothermal sink runs, 
the temperature potential was much greater than that for the 
low mass flow rate runs, thus the Grashof number was expected 
to be more significant in these processes. Since the Grashof 
number is greatly affected by a pressure increase as dictated 
by the p 2 term, it was felt that the Nusselt number would in- 
crease throughout a run. Once again these assumptions seem 
to be supported by the data, see figures 11, 12, 13 and 14 as 
well as 19, 20 and 21. 

As to the actual magnitude of the convective heat transfer 
coefficient as indicated by NTU, the data was hard to corre- 
late. As seemed consistent with the conclusions drawn above, 
the values of the Nusselt and Grashof numbers for the low mass 
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flow rate runs were lower than those for the high mass flow 
rate and isothermal sink cases. However, with the exception 
of the set of runs at a very low mass flow rate (perhaps cor- 
responding to the low AT discussion above) , the isothermal 
sink values were above those for the high mass flow rate runs. 
One possible explanation for this was the fact that the iso- 
thermal sink runs were recorded in 1 second intervals for run 
times of 7 to 9 seconds. The high mass flow rate runs , on 
the other hand, were only of 1 second duration with the points 
measured at 1/10 second intervals. Therefore it might be that 
for runs of larger duration, the mixing motion caused by the 
impinging jet could spread more thoroughly throughout the re- 
ceiver. These types of forced convection phenomena could not 
be adequately described on the basis of the results of these 
tests. In any case, the values calculated in these tests are 
at most 1 order of magnitude above those predicted by free 
convection theories (Figs. 11-21). Therefore if no better 
estimate of the heat transfer coefficient can be obtained, 
a value based on these theories could be used for processes 
of the nature discussed here. 

The preceding discussion was based mainly on physical 
intuition and is only crudely substantiated by the few experi- 
mental runs performed in this undertaking. As will be 
repeated in the recommendations, a much more extensive and 
controlled study is needed in this area to draw any conclu- 
sive relationships between the convective heat transfer co- 
efficient and the many variables in a system being charged. 
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VIII. CONCLUSIONS 



The general conclusions to be drawn from this investiga- 
tion are summarized as follows. 

(a) Excellent quantitative agreement between the 
experimental values and those predicted by Reynolds 
[Ref. 1] has been obtained for the thermodvnamic 
state of a gas in a receiver being charged. 

(b) The parameter NTU prescribed by Reynolds is a good 
indicator for describing the extent of heat transfer 
in a receiver being charged. 

(c) The technique developed in this paper for experi- 
mentally determining the convective heat transfer 
coefficient may be confidently used in conjunction 
with Reynolds' closed form solutions if the tem- 
perature potential between the receiver walls and 
the charging gas is of an order of 10°R or larger. 

(d) The convective heat transfer coefficient is related 
to the initial mass flow rate of the entering gas 
for all but the extremely slow charcring rates 

(Wq < .003 lbm/sec.). 

(e) The heat transfer coefficient for charging of the 
nature described in this paper can be estimated 
within one order of magnitude by free convection 
theories. The heat transfer coefficient thus 
obtained will be equal to or less than the actual 
value. The effects of errors in this estimate 
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will be small provided that the NTU is sufficiently 
small (high injection rates) . 

IX. RECOMMENDATIONS 

(1) It is recommended that an extensive and controlled study 
of the effects of the various system parameters on the 
convective heat transfer coefficient during charcring be 
made with emphasis placed on the development of an 
empirical formula for use in this area. 

(2) A possible extension of the technique for experimentally 
determining the convective heat transfer coefficient 
developed in this paper could include the use of analocr 
computer methods to directly analyze the pressure and 
temperature signals from probes in the vessel. 

(3) In any further investigation of the type discussed in 
this paper, the following suggestions are made. A more 
exact method of measuring the mass flow rate should be 
developed. Care should also be taken to eliminate any 
possible transient effects in establishing flows used in 
experiments of this nature. It is also recommended that 
hot gases be used in order to obtain substantial temper- 
ature potentials between the receiver walls and the gas. 
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APPENDIX A 



Nomenclature of Reference 1 

English Letter Symbols 
A Area, ft 2 

C Thermal Capacitance of receiver shell, Btu/°R 

c 

Cp Specific heat at constant pressure, Btu/(lb°R) 

c v Specific heat at constant volume, Btu/(lb°R) 

h Specific enthalpy , Btu/lb 

h Unit heat transfer convective conductance, 

Btu/(hr ft 2o R) 

k Thermal conductivity, Btu/(hr ft 2o R/ft) 

M Mass, lb 

P Pressure, lb/ft 2 

q Heat transfer rate, Btu/hr 

T Absolute temperature, °R 

U Total internal energy , Btu 

u Specific internal energy, Btu/lb 

V Volume, ft 3 

w Mass flow rate , lb /hr 

Greek Letter Symbols 
6 Time , hr 

y Viscosity, lb/ (hr ft) 

p Density, lb/ft 3 

A Denotes a difference 
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Nondimensional Grouping 



c* 


V (M 0 


c ) 

V 




k 


c /c 

p V 






M* 


M/M q 






NTU 


1/O^c 


W» ) 
V 0 


= (hA) . /c w„ 
v 0 


NTU 

CO 


V<R_c. 


W„) 

v 0 


= (hA) /c w n 
»' v 0 


p* 


P/P 0 






T * 


T/T 0 






T* 


T l /T 0 






T * 
c 


V T 0 






T * 

00 


T </ T o 






w* 


w/” 0 






Subscripts 






0 


Refers 


to 


initial conditions 


1 


Refers 


to 


inlet state 


c 


Refers 


to 


capacitance 


i 


Refers 


to 


inside receiver 


00 


Refers 


to 


environmental conditions 
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APPENDIX B 



Derivation of General Differential Equation and Closed Form 
Solutions for Charging 
Assumptions and Idealizations 

(1) The walls of the receiver are at a uniform temperature 
throughout . 

(2) The heat transfer resistances both inside and outside 
are constant and uniform throughout their respective 
surface areas. 

(3) Perfect mixing of the injected gas and the gas in the 
receiver. 

(4) The working fluid is a thermally and calorificallv 
perfect gas. 

(5) The specific heat of the capacitance is a constant. 

(6) The kinetic energy of the gas in the receiver may be 
neglected. 

General Differential Equation for Charging 
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An energy balance on the receiver, using h ^ as the stagnation 
enthalpy of the inlet gas gives , 



, dM _ dU 

h l de “ d0 q i ‘ 

An energy balance on the receiver walls yields, 

dU 



qi = 



de 



+ q 

” c 



The heat transfer rates can be defined in terms of the temper- 
ature potential and the heat transfer resistance as follows: 



T-T T -T 

c _ c °° 

q i ” R. ' q °° — R 

l 00 

The change in the total internal energy of the gas in the 
receiver may be written in terms of the mass and the specific 
internal energy as, 



dU ..du dM 

de " M de + u de * 

The rate of change of energy in the capacitance may be expres- 
sed as , 

dU d (C T ) dT 

c _ c c _ c 

de ” de c c de 

Thus, the energy balance on the gas may be written as, 

T-T 

7 _ „„du dM , c 

h i de “ M de u de + r. 

i 



or 



d (c T) T-T 

, 7 \ dM _ . . v ' c 

(h 1 u) dQ M d0 R . 



and thus , 



T-T 

, m m\dM ^dT , c 

(c T.-c T)-^ = C M-r^- + — 

pi v de v de R- 

^ l 
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Note here that is the stagnation temperature of the inlet 
gas. Therefore, dividing by gives, 

(kT _ T )^i - 

(KT 1 T; de M d0 R. 



In a similar fashion the energy balance on the capacitance 
can be written as , 

T-T dT T -T 

£ = c + c 00 



'c de 



then dividing by C c , 



dT T-T T -T 
c c c °° 



d6 



R C R C 
ic °° c 



There fore , 



= JR + — + i_ ) 

d6 R.C R C C y R. R 

1 C oo c Cl 00 



The instantaneous mass flow rate, w, may be defined as 

dM 



w 



d6 



and we can write. 



as well as, 



dT dT dM dT 

d0 dM d0 W dM 



dT dT ... dT 
c _ c dM _ c 

de " dM d0 " W dM 



The equation for the energy of the gas may then be written as, 

<kT -T)w = WMJ§ + ^ 

1 V 



or 



T-T 

WM ® + iC^T + * 0 • 

1 V 



( 1 ) 
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In a similar manner the energy equation of the caoacitance 
can be written as 



w- 



dT 
c 

dM 



R i C c 



T T , 1 

" _£(i_ + L_) 

R C C V R. R ' 

oo c Cl 00 



or 



dT T , , 

c , c / 1 , 1 > T 

dM C V R. R ' R.C 

Cl 00 1C 



'00 = 0 . 



R C 

oo c 



( 2 ) 



Expressing equation (1) in dimensionless form by dividing 
through by and noting that 



M^— = — — — — = —5 

dM M ft , ,M * divF 3 

0 dfjj-O 



we obtain, 



T * — T* 

,„*dT* , L c 

w M*__x- + 

dM* R. c 
l v 



+ (T*-kT*)w = 0 . 



Now dividing by w Q and noting the definition of NTU, we ob- 
tain the dimensionless form, 

rlT* 

w*M*5^r + (T*-T* ) NTU + (T*-kT*)w* = 0 . (la) 

dM* c 1 

In the same way, equation (2) can be non-dimensionalized by 
dividing through by TqWq/Mq so that, 



w 



. dT S , T ? M o.i , 1 , , T * M 0 T i M 0 

dM* C c w 0 ' Rl R J Ri w 0 C c - R;o C c w 0 



= 0 



and noting that NTU/C* is equal to ’^o //R i W 0 C c t ^ le 



is obtained, 



dT* 

c 

dM*" 



w* 



+ T* (NTU+NTU ) - 



NTU, 

C* 



NTU 



•T* - 



Cj— co 



■T* = 0 . 



(2a) 



These two dimensionless equations can be combined to yield a 
single equation giving T* as a function of M* by differentiation 
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of equation (la) with respect to M* and suitable arithmetic 
manipulations represented below: 

.dT* 



+ (T*-T*)NTU + (T*-kT*)w* = 0 



(la) 



Differentiation with respect to M* yields. 



dw* dT* 

dM*" M + dM* 



) + w*( 



dT* 

dM* 



+ M* 



d 2 T* 
dM* 2 



) + ( 



dT* 



dT* 

c 



dM* dM 3 



) NTU 



+ + ( S£ )t * - kT i ( ^ - 0 

rearranging) the following is obtained. 



dT* 

c 

dM* 



*,~dT* . * d 2 T* 

W 2 dM* + m dM* 2 dT* 
NTU + dM 7 



dw* /m*. w*dT*\ 
— r* (t*+ m *^tt) 



dM’ 



dM’ 



NTU 



w*kT*^^— 
W Ki ldM* 

NTU 



Substituting this expression for dT*/dM* in equation (2a) , 



* 2 *d 2 T* 

w m *m*r 

NTU 



2w* 2 dl P* 

W dM* , „*dT* 
NTU + W dM* + 



,dw’ 



.dT*. 



w 'dS* (T * +M *atV ) 

NTU 



W*kT*^*- T* NTU 

- + — — (NTU+NTU ) - T* - — T* = 0 

NTU C * U C* C* “ 

Multiplying through by NTU and combining gives , 

d 2 T * dw* dT* 

w * M *|pr + w*(2w*+NTU+M*§^ r )^ r 



+ ( W *^l _ NTU 2 * _ 
+ ' w dM* C* 1 



NTU NTU d * 

-T* - w*kT*- 



C* 

^0 



ldM* 



( 3 ) 



T * 

+ ^-(NTU+NTUjNTU = 0 
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Now solving (la) for T* , 



dT* 

w*M*^r + T*NTU+T*w*-kT*w* 
T * = dM* 1_ 

c NTU 

Therefore the term, 

T* 

(NTU+NTU^) NTU 



in the above expression can be written as, 



dT * NTU+NTU, 
W * M *dM* *■ C* ^ 



T*NTU 

C 0 



(NTU+NTU,) 



w *t* kT l W * 

+ (NTU+NTU ) - — — (NTU+NTU ) 



C* 



'0 



Substituting this in equation (3) we arrive at the general 
differential equation for charging. 



w*M* 



d 2 T* 
dM* 2 



+ [w 



rdw* 

dM* 



d * ntu+ntu to 

w* t2w*+NTU+M*^- + ( ^ — ) M* ] 

, NTU NTU, 

+ ^r(NTU+NTU )w* + 7T* — ] T* 

c o 00 C 0 



dT* 

dM 7 



- VT*w*— 
Kr i W dM* 



NTU+NTU. 
C* 



-w*kT* - 



NTU NTU 

O 

c 0 



T* = 

00 



Closed Form Solution for Adiabatic Charging 



/ /ll dU - d(Mu) 

/z^dM 

/77/77/7T/ 
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Since there are no heat transfer rates to consider in an 
adiabatic solution it is not necessary to specify the mass 
flow rate. An energy balance on the gas in the receiver, 
using stagnation functions, gives 

Ti^dM = d(mu) = Mdu + udM 



or 



(Ti^-u)dM = Mdu 



Then one can write, 

(c T, -c T)dM = Me dT . 
p 1 v v 

Dividing by c^ gives , 

(k^-T) dM = Mdt . 



In order to get this equation in dimensionless form, we 
divide through by TqMq, thus obtaining 

(kT*-T*)dM* = M*dT* . 



Now separating variables and integrating this equation from 
the initial conditions (T = Tq, M = Mq) to a state at a later 
time (T = T, M = M) , 




This gives , 

-logtk^-T*) + log(kT*-l) = logM* 



Rearranging, this equation reduces to 



(kT*-l) 

M* 



kT* - T* 
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Solving for T* , the final result is obtained. 



T* = kT* 



(kT*-l) 



Thus having derived a relation for T* as a function of M* , 
the perfect gas relation P* = M*T* is used to obtain, 

P* = kT*(M*-l) + 1 

Closed Form Solution for Isothermal Charging 




Since T = constant = T Q = T c , then 




Therefore the isothermal pressure mass relationship can be 
obtained directly from the perfect gas law, and is 

P* = M* 
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Closed Form Solutions for Charging at Constant Mass Flow with 
Heat Transfer to an Isothermal Sink 




An energy balance on the gas using stagnation functions gives, 



dM 

J ld0 



= q,- + 



d (Mu) 

de 



dM . du. 



" + Ujg- + 



By definition of the heat transfer resistance R^ , 

T-T 



q i = — 



and using w = dM/d0 , then 

T-T 



C P T 1 W 0 = — 



£ + M 3F + w 0 c v T 



or 



T-T 



c T, w A = 



c . dT , m 

— + Me -T7T + w A c T 
vd6 0 v 



pl w 0‘ R 

Using the following relation, 

dT dT dM dT 

d9 dM d0 W 0dM 

and dividing by Wq and rearranging the exoression becomes, 

T-T 

dT c 

-c T. + c M-^tt + c T + = 0 

pi v dM v R.w A 

c l 0 



71 



Now dividing by to obtain 



T-T 

dT g 

m— + T - kT. + - — — 
dM 1 R.w rt c 



= 0 



i 0 v 

Note that all functions refer to stagnation conditions. In 
order to get this equation in dimensionless form it is divided 
by Tq so that, 

_l fp* _ VT 1 * 4- t'V*— T*'l 1 

T„ dM + 1 U c'r.w.c 

0 l 0 v 



Since 



M dT 
T q dM 



d(|-) 

M 0 



M o d<Sr> 

0 



dT* 1 

M*^tW- and NTU = 

dM* R . w 



0 c v 



this expression can be written as, 
dT* 

M*^-=- + (T*-kT*) + (T*-T* ) NTU = 0 . 

dM* 1 c 

In order to obtain the closed form solution the variables 
are separated and the expression integrated from the initial 
conditions (T = Tq , M = Mq) to a state at a later time (T = T, 
M = M) . 



( 




dT* 

[T* ( 1+NTU) -kT*-NTUT* ] 
i. c 




M 

M 0 dM* 




This gives 



r r +^ ' T U ) log [-kT*-NTUT*+ (1+NTU) T*] 

- (I+ - ^ TU ) - log t-kT*-NTUT*+ (1+NTU) ] = -logM* 
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or 

-kT* - NTUT* + (1+NTU) T* = M * _tl+NTU ^ [-kT, -NTUT*+ ( 1+NTU) ] . 

-L C 1 C 

Solving this expression for T* as a function of M* the follow- 
ing is obtained: 



T* = 



kT?+NTUT*- (kT*+NTUT*-l-NTU) M* 
1 c 1 c 



- (1+NTU) 



(1+NTU) 



Using the perfect gas law. 



P* = 



(kT*+NTUT*)M*- (kT*-l-NTU+NTUT*)M* 

1+NTU 



-NTU 



Closed Form Solution for Charging at Constant Mass Flow with 
Inside Resistance Negligible 




An energy balance on the gas and the capacitance combined 
yields , 



rlT 

U dM _ d (Mu) , n < 

ide de ^cde 



+ q 



By definition of the heat transfer resistance between the 
capacitance and the surrounding medium, R , 



Sco = 



T -T T-T 

Q 00 00 



R 



R 
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Now since Wg = dM/d0 the expression 



dT _ dT dM dT 

d0 dM d0 w 0dM 

can be written. Using this expression and the perfect aas 
relationships in the original energy balance, 

T-T 



m dM , dM _ , dT dT 

— c T. -577 + . - c T + Me + C -377 + — =— 
p ld0 d0 v vd0 cd0 R 

c Cx 



= 0 



Dividing through by c^dM/d© or c^Wg and noting that dT c /d0 
= dT/d0 , 

C T-T 

r* /-J r P no 

= 0 . 



C T-T 

c 00 

-kT. + T + (M+— )=±- + — = 

1 c dM C R w„ 

V V 00 0 



In order to express this equation in terms of non-dimensional 
temperature, T* , it is divided by Tg to obtain, 
dT* dT* 

M*^T* + C*2i^ _ kT* + T* (T*-T*) NTU = 0 . 

dM* OdM* 1 co 00 

The variables are now separated and the expression integrated 

from the initial conditions (T = T« , M = M n ) to a state at a 



later time (T = T, M = M) . 



T_ 

T, 



M_ 

M, 



dT* 

T* ( 1+NTU ) -kT*-NTU T* 



dM’ 



Mg ^ +l 0 

mT 



The solution to the above gives, 



(M*+Cg ) 



log [T* (1+NTUj -kT*-NTUT*] = - (1+NTUj lo( ? JX+C* 7 
+ log [-kT*-NTU eo T*+ (1+NTUj ] 
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or 



T* ( 1+NTU ) - kT* - NTU T* 

00 00 00 



(M*+Cq) 



TI+C*T 



- (1+NTUJ 

[-kT*-NTU x) T£+ (1+NTU) ] 



Solving this expression for T* gives 



/C*+l 



yl+NTU 



T * “ 



(1+NTU -kT*-NTU T ) C*+M*l + kT* + NTU T* 

CO J_ CO GO \ 0 / 1 0000 



1+NTU 

00 

This expression may be combined with the perfect gas law in 
order to obtain an expression for P* as a function of M*. 

1+NTU 



>* - 



M* 



c 3 +1 

(1+NTU -kT*-NTU T ) C*+M*l + kT*+NTU T* 

OO 1 oo 00 | 0 / 1 “ “ 

1+NTU 
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Table 1 



Physical Dimensions of Experimental Apparatus 

1. Test Receiver 

Volume: 1.03 cubic feet 

Internal Area: 5.61 square feet 

Nominal Internal Diameter: 10.19 inches 

Nominal Outside Diameter: 10.75 inches 

Nominal Internal Height: 21.75 inches 

Mass of the Receiver: 38.5 lbms . 

Material: 6061-T6 Aluminum 

Specific Heat: .226 Btu/lbm. °F 

Thermal Conductivity: 126 Btu/hr. ft. °F 

2. Flow Metering Devices 

3/8 in. knife-edge orifice - discharge coeff. = .620 
1/8 in. nozzle - discharge coeff. = .858 
1/32 in. nozzle - discharge coeff. = .856 
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Table 2 





S uirana ry 


of High Flow 


Rate Runs 






C* = 657 


k = 1.4 


T* = 1.00 




Run 


Total Run 


Mass Flow 


Average 


Average % P* 


No. 


Time 


Rate 


NTU 


Deviation 




(sec) 


(Ibm/sec) 






1 


1 


.327 


.123 


7.93 


2 


1 


.325 


.092 


8.11 


3 


1 


.325 


.104 


5.59 


4 


1 


. 310 


.101 


7.47 


5 


1 


.309 


.089 


6.93 


6 


1 


. 309 


.085 


6.11 


7 


1 


. 294 


.121 


10.06 


8 


1 


.293 


. 084 


5.58 


9 


1 


.293 


.084 


4.92 


10 


1 


.278 


. 102 


6.93 


11 


1 


.277 


. 105 


7.35 


12 


1 


.277 


.115 


9.63 
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Table 3 



Summary of Low Flow Rate Runs 



Run 

No. 


Total Run 
Time 
(sec) 


Mass Flow 
Rate 

(lbm/sec) 


Average 

NTU 


Average % P* 
Devi ation + 


1 


105 


. 0031 


3.99 


-2.61 


2 


105 


.0031 


2.86 


-3.41 


3 


105 


.0031 


2.90 


-3.18 


4 


105 


.0029 


3.69 


-2.84 


5 


105 


. 0029 


3.01 


-3.11 


6 


105 


. 0029 


4.78 


- 2.71 


7 


105 


.0028 


1.86 


-3.68 


8 


105 


.0028 


3.49 


-2.93 


9 


105 


.0028 


2.18 


-3.39 


10 


105 


.0024 


2.92 


- 2.76 


11 


105 


.0025 


3.39 


-2.65 


12 


105 


.0025 


5.54 


-2.22 


13 


105 


.0022 


5.34 


-2.50 


14 


105 


.0022 


7.78 


-1.97 


15 


105 


.0022 


5.61 


-2.09 


16 


105 


.0018 


4.52 


-2.02 


17 


105 


.0018 


7.80 


-1.67 


18 


105 


.0019 


4.43 


-2.01 


19 


105 


.0015 


5.17 


-1.70 


20 


105 


.0015 


8.00 


-1.44 


21 


105 


.0015 


7.91 


-1.27 


22 


105 


.0012 


4.89 


-1.54 


23 


105 


.0012 


6.02 


-1.33 


24 


105 


.0012 


3.83 


-1.70 


25 


105 


. 0009 


4.32 


-1.43 


26 


105 


.0009 


6.01 


-1.12 


27 


105 


.0009 


6.21 


-1.07 


28 


105 


.0007 


4.58 


-1.49 


29 


105 


. 0007 


3.77 


-1.36 


30 


105 


. 0007 


3.47 


-1.45 


31 


90 


.0006 


4.45 


-1.07 


32 


90 


. 0006 


13.88 


-0.46 


33 


90 


.0006 


11.19 


-0.52 


34 


75 


.0004 


14.49 


-0.54 


35 


75 


.0004 


14.62 


-0.54 


36 


75 


.0004 


12.39 


-0.51 



+ A negative deviation indicates that the P* values predicted 
by Reynolds were lower than those actually observed. 
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Table 4 





S urrana rv 


of Isothermal 


Sink Runs 




Run 


Total Run 


Mass Flow 


Average 


Average % P 


No. 


Time 


Rate 


NTU 


Deviation 




(sec) 


(lbm/sec) 






1 


150 


.0031 


2.69 


-10.39 


2 


150 


. 0031 


2.24 


- 8.84 


3 


150 


.0031 


2.17 


-11.08 


4 


150 


.0028 


2.07 


-11.59 


5 


150 


.0028 


2.41 


-10.35 


6 


150 


.0028 


1.97 


-11.47 


7 


150 


. 0025 


3.02 


- 9.10 


8 


150 


.0025 


3.30 


- 8.58 


9 


150 


.0026 


2.94 


- 8.94 


10 


9 


.0488 


1.18 


-13. 32 


11 


9 


.0478 


1.41 


-13.40 


12 


9 


.0486 


1.38 


-13.22 


13 


9 


.0420 


1.39 


-12.67 


14 


9 


. 0422 


1.28 


-13.37 


15 


9 


.0422 


1.15 


-14.34 


16 


9 


.0350 


1.41 


-12.49 


17 


9 


.0348 


1.03 


-13.68 


18 


9 


.0349 


.94 


-13.57 


19 


9 


.0274 


1.73 


-10.78 


20 


9 


.0272 


1.13 


-12.36 


21 


9 


.0272 


1.20 


-12.47 


22 


9 


.0194 


1.85 


- 9.92 


23 


8 


.0192 


1.09 


-11.14 


24 


8 


.0190 


1.05 


-11.14 


25 


7 


.0111 


1.29 


- 7.68 


26 


7 


.0111 


1.20 


- 7.85 


27 


7 


.0114 


1.15 


- 7.98 



+ A negative deviation indicates that the P* values oredicted 
by Reynolds were lower than those actually observed. 
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+ * 



SCHEMATIC OF TEST RECEIVER 
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2-pressure transducer 1-thermocoupl e patch panel 

TEST SYSTEM SCHEMATIC 
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Photograph of Test Receiver 
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Photograph of Test System 
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2 = x 

HIGH MASS FLOW RUNS 1-3 3 = * 

( 1 ) laminar free convection line 

(2) turbulent free convection line 

Pig. 11 
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(1) laminar free convection line 

(2) turbulent free convection line 

Fig. 12 
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8 = x 

HIGH MASS PLOW RUNS 7-9 9 = * 

(1) laminar free convection line 

(2) turbulent free convection line 

Pig. 13 
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GRPR 10 = + 

1 1 = x 

HIGH MASS PLOW RUNS 10-12 1 2 = * 

(1 ) laminar free convection line 

(2) turbulent free convection line 

Pig. 14 
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2 = x 

LOtf MASS PLOW RUNS 1-3 3 = * 

1 ) laminar free convection line 

2) turbulent free convection line 

Pig. 15 
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(1) laminar free convection line 

(2) turbulent free convection line 

Fig. 16 
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(1) laninar free convection line 

(2) turbulent free convection line 



Pig. 17 
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1 4 = x 

LOW MASS FLOW RUNS 13-15 1 5 = # 

(1) laminar free convection line 

(2) turbulent free convection line 

Fig. 18 
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( 2 ) 



,0 8 


10 9 


101° 




3RPR 


4 = + 






5 = x 




ISOTHERMAL SINK RUNS 4-6 


6 = * 



(1) laminar free convection line 

(2) turbulent free convection line 

Pig. 19 
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ORPR 10 = + 

1 1 = x 

ISOTHERMAL SINK RUNS 10-12 1 2 s * 

(1) laminar free convection line 

(2) turbulent free convection line 

Pig. 20 
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GRPR 16 = + 

17 = x 

ISOTHERMAL SINK RUNS 16-18 1 8 = * 



1 ) laminar free convection line 

2) turbulent free convection line 

Pig. 21 
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